Doppler cooling of an optically dense cloud of trapped atoms 
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We have studied a general technique for laser cooling a cloud of polarized trapped atoms down 
to the Doppler temperature. A one-dimensional optical molasses using polarized light cools the 
axial motional degree of freedom of the atoms in the trap. Cooling of the radial degrees of freedom 
can be modelled by reabsorption of scattered photons in the optically dense cloud. We present 
experimental results for a cloud of chromium atoms in a magnetic trap. A simple model based 
on rate equations shows quantitative agreement with the experimental results. This scheme allows 
us to readily prepare a dense cloud of atoms in a magnetic trap with ideal starting conditions for 
evaporative cooling. 



O 
O 

(N 

< 
o 

(N 



> 
On 
(N ■ 

oo : 
o . 

(N ■ 

o ■ 

^— > ■ 
G ■ 
cd 

3 : 



PACS numbers: 020.4180, 020.7010, 140.3320 



I. INTRODUCTION 

Since the advent of laser cooling of neutral atoms more 
than 20 years ago, this method has been extensively stud- 
ied |lj and widely used for a variety of atoms and appli- 
cations. Ultimately, laser cooling followed by evaporative 
cooling Q| allowed the creation of a Bose-Einstein con- 
densate (BEC) of neutral atoms 

In almost all BEC experiments atoms are caught in a 
magneto-optical trap. After sub-Doppler molasses cool- 
ing Q and polarization of the sample by optical pump- 
ing, the atoms are typically transferred into a magnetic 
trap H, H> |§L or , as has been recently demonstrated, into 
an optical trap. The consecutive evaporative cooling 
process requires both, a high initial density to provide a 
reasonably large elastic collision rate, and a large num- 
ber of atoms, since most atoms are removed from the 
trap during evaporation. The figure of merit for efficient 
evaporative cooling is the initial phase space density of 
the atomic cloud. Therefore, all laser cooling and polar- 
ization steps have to be optimized carefully, to keep the 
temperature of the atoms low and the number of atoms 
high, thus maximizing the phase space density in the final 
trap. 

In this paper, we present a robust Doppler cooling 
scheme in an external trap which is applicable to most 
laser-coolable atomic species. Our scheme is particularly 
well suited for optically dense samples and reduces the 
figure of merit for evaporative cooling to the number of 
atoms transferred into the trap regardless of tempera- 
ture. 

Free-space one-dimensional Doppler cooling is usually 
performed in a standing wave light field created by two 
counterpropagating laser beams with a frequency below 
the atomic transition frequency. Cooling is based on pref- 
erential absorption of photons from the laser beam op- 
posing the direction of motion of the atom. Subsequent 
spontaneous emission is centrally symmetric and does not 
change the mean momentum of the atom. Doppler cool- 
ing of polarized atoms in a magnetic trap has been pro- 
posed for atomic hydrogen Q and experimentally real- 
ized for sodium ||, hydrogen [jl0| and lithium JO]. In 



Ref. H one-dimensional Doppler cooling of an optically 
thin cloud of sodium atoms in a magnetic trap was per- 
formed. The atoms were treated as a two level system 
since a very high magnetic offset field (Bq ~ 1500 G) 
spectroscopically resolved the Zeeman substates. Cool- 
ing of the motional degrees of freedom orthogonal to the 
laser beams was provided by anharmonic mixing in the 
trapping potential. Due to a long mixing time compared 
to the cooling time, the achieved temperature was ten 
times the Doppler temperature for sodium. In Ref. [h] 
single-beam pulsed Doppler cooling was performed on a 
dense sample of magnetically trapped hydrogen atoms. 
Thermalization of the cloud was accomplished by elastic 
collisions between the atoms. Cooling was limited to ap- 
prox. five times the Doppler temperature due to limited 
laser power and additional heating rates. 

The main advantage of the cooling technique presented 
here is that it combines three-dimensional temperatures 
close to the Doppler-limit with the cooling of dense sam- 
ples. In our experiment, one-dimensional optical mo- 
lasses cools the cloud of atoms in axial direction down 
to the Doppler limit ji]. Cooling in the radial directions, 
orthogonal to the laser beams, can be explained by re- 
absorption of spontaneously emitted photons by the op- 
tically dense cloud. To our knowledge, reabsorption has 
so far been treated only as a density-limiting mechanism 
which comes along with heating of the atomic sample 
p"2| , fuH |l4| , p^f . Here, we focus on the cooling aspects of 
scattered and reabsorbed photons in a trapped polarized 
atomic sample. 

We show that the atoms in the magnetic trap remain 
polarized during the cooling process and practically no 
atoms are lost, provided the magnetic substructure is 
spectrally resolved. We are able to cool chromium atoms 
in a magnetic trap from ps 1 mK to a mean temperature 
of 240 fiK, corresponding to an increase in phase space 
density by a factor of 80. We have studied the dynamics 
of the cooling process as well as the steady state temper- 
ature for various cooling parameters. The experimental 
results can be explained with a simple model based on 
rate equations for the temperature. 

The paper is organized as follows. In Sec. 2, we present 
a theoretical model for Doppler cooling in a magnetic 
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trap in the presence of reabsorption. The setup and ex- 
perimental techniques for our cooling experiments with 
chromium are described in Sec. 3. Results on the dy- 
namics and steady state properties of the cooling process 
are presented in Sec. 4. We conclude with a discussion 
of our results in Sec. 5. 



II. THEORY 

A. Rate Equations 

In this section, we present a simple model based on 
rate equations. Although we assume cooling in a mag- 
netic trap here, the scheme is universal to all kinds of 
traps, provided cooling is compatible with the trap and 
the quantization axis of the atoms is independent of the 
position of the atoms in the trap. 

Consider a cloud of spin-polarized atoms with mass 
m and magnetic moment /i confined in an axially sym- 
metric harmonic magnetic trap with trapping potential 
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FIG. 1: Part of the Zeeman substructure of 52 Cr as an exam- 
ple of a J — > J + 1 transition. Shown is the transition from 
the electronic ground state 7 S3 to 7 P4 with a wavelength of 
Aatom = 425.55 nm, a saturation intensity I s = 85.2 W/m 2 
and a linewidth of V = 2tt x 5.02 MHz. The numbers next to 
the transitions are the squares of the Clebsch-Gordan coeffi- 
cients. 



Zeeman-splitting due to the finite offset field B addi- 



B are the magnetic field curvatures in radial (x, tionaUy reduces ^ scattering rates of unwanted 



y) and axial (z) direction, respectively, and Bq is a ho 
mogeneous offset field along the axial direction. We as 
sume that the corresponding trap frequencies u) x .y,z = 

^ Bx - V - Z are much smaller than the Larmor precession 

frequency at the center of the trap lol — ensuring 
that the atoms stay spin polarized and no spin flip losses 
occur [ ^6[ . 

Doppler cooling is performed on an electric dipole tran- 
sition that couples a long-lived ground state with total 
angular momentum J to a short-lived excited state with 
angular momentum J , such that J = J + 1 jllj . We 
also assume the corresponding Lande factors gj and g y 
to have the same sign, which we choose to be positive for 
our example. The atoms are polarized in the "weak field 
seeking" state mj — J and have a magnetic moment of 
(i = mjgjfiB, where /is is Bohr's magneton. 

Cooling is carried out in a a + /er + -standing wave cre- 
ated by two laser beams propagating along the axial 
direction with normalized intensity / = /laser /-^sat and 
detuning A = (uvaser — ^atom)/r from the unperturbed 
atomic transition frequency. / sa t and T are the saturation 
intensity and the linewidth of the transition, respectively. 

In general, the quality of the polarization of the cool- 
ing light is important. For perfectly polarized er + -light 
propagating along the quantization axis of the atoms in 
a homogeneous magnetic field, no depolarizing transi- 
tions mj — ► mjjinj — 1 can occur. This is also true 
for excitation by scattered photons since the polariza- 
tion with respect to the quantization axis of the incident 
light is preserved in the scattering process if the prop- 
agation direction coincides with the orientation of the 
magnetic field. For imperfect polarization, the cooling 
light can also drive depolarizing transitions which cause 
loss of atoms. In the configuration described above, the 



depolarizing transitions with respect to the rate i? po i of 
the cooling transition mj — ► mj + f : 
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where A po i and Aa ep are the effective detunings between 
the laser and the atomic transition frequency in a mag- 
netic field and Cf is the square of the Clebsch-Gordan 
coefficient for the transition from state m g = i to m e = j. 
In Fig. [I] this effect is illustrated for the J = 3 — > J = 4 
transition that connects the ground state 7 S3 to the ex- 
cited state 7 P4 in 52 Cr. A full treatment of optical pump- 
ing in an arbitrary light field taking into account all ex- 
cited and ground state levels gives the steady state pop- 
ulation of each ground state for a given magnetic field 
|l8| . Assuming a polarization quality of 1:10, a detuning 
of A = —0.5 for the polarizing transition and allowing for 
an atom loss of less than 1% during the cooling process, 
we find for the case of chromium a necessary offset field 
of Bq > 5 G and for rubidium Bq > 15 G. 

Following the common treatment of free-space Doppler 
cooling p], El, rate equations for the temporal evolution 
of the temperature for the axial (z) and radial (y) degrees 
of freedom are readily obtained: 
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where Er = *Ujr is the recoil energy of the cooling transi- 
tion with wavenumber k — 2tt/X and k B is the Boltzmann 
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constant. Usually, the heating and cooling time constants 
are derived by considering the random walk in momen- 
tum space due to spontaneous emission and the net cool- 
ing effect of scattered photons from the two counterprop- 
agating laser beams \%. In an optically dense sample 
reabsorption of photons has to be taken into account. 
Previously, reabsorption has mainly been discussed as 
a reason for density limitations and density-dependent 
heating in magneto-optical traps ]l2| , fi"3[ |l4| and in the 
context of sub-recoil cooling p5| . In our model, we em- 
phasize cooling effects of reabsorbed photons which can 
be substantial for polarized atoms confined in an external 
trapping potential. 

The photon energy after a scattering event in the low 
intensity limit |ll| is shifted by at most AEr + 2hkv <C 
7i|A P oi|r, where v is the velocity of an atom at Doppler 
temperature. Therefore scattered photons have essen- 
tially the same detuning and polarization as the cooling 
laser beams. For our calculations, we use the axial direc- 
tion parallel to the laser beam as our quantization axis. 
The contribution of the reabsorbed photons to cooling of 
the radial and axial directions depends on the projection 
of their wave vector k onto these directions. The effec- 
tive intensity I z s y of the cooling photons is proportional 
to the incident laser light intensity 21 with a proportion- 
ality factor K ZiV which primarily depends on the optical 
density OD zy in these directions: 
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If we consider a dipole radiation and absorption pattern 
for CT-polarized light and take into account the 1/e 2 size 
a{9) of a Gaussian shaped cloud in the direction given 
by 6, we arrive at the following expressions: 
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k z = k J d6 [l + cos 2 ff) 2 a {6) sin 6\ cos 6\ (7) 
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kq = tiq(J\ — ^-j= is an angle-independent constant, with 
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peak density no and resonant absorption cross section 
a\ = 6ttX 2 which we assume to be uniform over the cloud. 
In a simple picture, k, ZiV can be interpreted as the effec- 
tive number of photons that contribute to cooling in the 
corresponding directions originating from a single photon 
absorbed from the laser light. The total number of reab- 
sorption/emission cycles for an incident photon is given 
by 



K = k J dO(l + cos 2 0) 2 <t{9) sint 
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In axial and radial direction the total cooling light inten- 
sity reads l\ ot = 21 + I v z = 21(1 + k z ) and J* ot = 1% = 
21 k v , respectively. 

For the derivation of the heating rates, we have to take 
into account directed reabsorption of K z y photons and 
spontaneous emission of (k + 1) photons for each photon 
absorbed from the cooling laser beams: 
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The factors 2/5 and 3/10 are statistical weights from 
the dipole radiation pattern for the different directions of 
spontaneous emission. In this model we have neglected 
additional loss and heating processes, like e.g. radiative 
escape ^(J and radiation trapping [Q. The influence of 
these additional effects will be discussed in the experi- 
mental part of this paper. 



B. Discussion of the Model 

1. Steady State 

For simplicity we assume the coefficients n z ^ y to be 
time independent. This is justified for the evaluation of 
the steady state parameters of the cooling process. The 
dynamics of the cooling process can be described only 
qualitatively within this model, since the optical density 
will change during cooling and therefore the number of 
reabsorbed photons increases. 

The steady state temperatures are readily derived from 
Eqns. § | to be 
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where To = jj^ is the Doppler temperature. These 
minimum temperatures are found at a detuning of A po i = 
— 5, in agreement with free space Doppler cooling [Q. 
It is worthwhile to mention that in this configuration 
in absence of reabsorption (i.e. n z ^ y = k = 0) the ID 
temperature in axial direction 



If the effects of reabsorbed photons are included the cool- 
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FIG. 2: Temperatures in axial (z) and radial (y) direction 



as a function of trap aspect ratio for N 
density of n — 5 x 10 10 cm -3 . 



10 atoms a peak 



is smaller than what is usually called the ID Doppler 
temperature. The reason is easily seen: in the ID model 
photons can be absorbed and emitted only along the axial 
direction; in a 3D model with ID cooling, photons can 
be emitted into radial degrees of freedom according to 
the dipol radiation pattern. These degrees of freedom 
take up photon recoil heating and therefore reduce the 
heating rate in the axial direction. In this situation the 
radial energy grows without limit. 

If we now include reabsorption, the situation changes. 
Whereas reabsorption only in axial direction (k 2 > 0) 
changes the cooling rate but not the steady state temper- 
ature, reabsorption in radial direction (n y > 0) provides a 
cooling effect, that establishes a radial steady state tem- 
perature (Eq. |I4| ). The magnitude of this temperature 
is determined by the ratio of heating to cooling photons. 

Experimentally, it is not possible to independently 
change and k. They depend via Eqns. [t] and 

H on the shape of the trapped cloud. In Fig. |^ we show 
an example, where we have integrated these equations for 
different aspect ratios a = o z j a v to calculate the steady 
state temperatures. The aspect ratio has no influence 
on the axial temperature, since cooling is dominated by 
photons scattered directly from the laser beams. In ra- 
dial direction we observe a minimum near a spherical 
cloud. For a cigar shaped cloud the optical density in ra- 
dial direction decreases and for a pancake-shaped cloud 
the solid angle decreases in which photons contributing 
to radial cooling can be scattered. 



2. Dynamics 

The transient evolution of the temperature in radial 
and axial direction can be approximated by the well 
known solution of Eqns. |^: 



where T® y are the initial temperatures of the cloud in 
axial and radial direction. Depending on the initial tem- 
peratures and the steady state temperatures Eqns. |l^, 
ful , cither Doppler cooling {T^° y <T° y ) or Doppler heat- 
ing {T^ y > T® y ) occurs. However, since the atoms are 
oscillating in a harmonic trap, cooling preferentially oc- 
curs close to the trap minimum where the velocity of the 
atoms is highest. Therefore the cooling rates can not be 
higher than approx. one quarter of the trap oscillation 
period t z ,y- 

According to Eqns. ^, |l(] the time constants rf° l , 
Ty°° l depend on the number of scattered cooling pho- 
tons. This can provide a test for the reabsorption model. 
If cooling in radial direction originated from elastic scat- 
tering pX| ] or anharmonic mixing by the trapping poten- 
tial ||, one would expect an intensity independent re- 
laxation towards an equilibrium temperature, provided 
T cooi ^ Ty°° l , while our model explicitly requires an in- 
tensity dependence of Ty°° l 

We have already mentioned in the previous section that 
the effective intensity of reabsorbed photons will change 
during the cooling process. We will now discuss the dy- 
namical effects on the relevant parameters K z . y . From 
Fig. |^ and Eq. [l^ we deduce a weak dependence of the 
axial steady state temperature on reabsorbed photons. 
In this regime cooling in axial direction is dominated by 
the real laser beams and we can assume T z to be unper- 
turbed from its minimum value T™ m . On the contrary, 
radial cooling strongly depends on the reabsorbed pho- 
tons, so we have to take into account how n v changes 
with time. As cooling in the z-direction starts, the cloud 

shrinks according to a z = A / p 3f ■ That will increase the 
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optical density in radial direction and thus the number 
of reabsorbed photons according to Eq. || and 

OD y oc oc Ky. 

O z Oy 

As a result, the radial steady state temperature is re- 
duced, since it depends on the effective number of re- 
absorbed photons (Eq. |l4|). The same argument holds 
for the size of the cloud in radial direction. During ra- 
dial cooling a y shrinks and we arrive at an even lower 
T^°. In actual experiments this nonlinear effect is less 
pronounced and will be neglected. 



III. EXPERIMENTAL SETUP 

In our experiment, we continuously load typically 
10 s chromium atoms in the low field seeking | j = 3, mj = 
3 > magnetic substate into a weak Ioffe-Pritchard trap 
pl| in the cloverleaf configuration p^j. Details of the 
continuous loading mechanism can be found in Refs. 
|23| , |24| . In the compressed cigar-shaped trap with a nearly 
harmonic trapping potential in all three dimensions, we 
have an axial and a radial curvature of B z =110 G/cm 2 
and B = 750 G/cm 2 , respectively. This corresponds 
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to trapping frequencies for chromium /i = 6/ib of w x = 
Lu y = 2n x 110 Hz and uj z = 2n x 42 Hz. We use a 
rather high offset field of Bq — 28 G to assure sufficient 
harmonicity in all directions. During compression the 
atoms are adiabatically heated to a temperature of 1 mK 
at a peak density of 0.4 — 1 x 10 10 cm' 3 . In this trap we 
perform one-dimensional Doppler cooling with a retrore- 
flected cr + -polarized beam along the axial direction of 
the trap. Laser light at 425.55 nm from a frequency dou- 
bled Ti:Sapphire laser drives the 7 S3 — > 7 P4 transition 
in 52 Cr (r = 2tt x 5.02 MHz, I sat = 85.2 W/m 2 , see Fig. 
|l|) . The laser frequency is tuned to 7 T above the unper- 
turbed transition frequency corresponding to an effective 
detuning of approx. —0.8 T for a nonmoving atom placed 
at the minimum of the magnetic field in the center of the 
trap. This value was optimized with respect to minimum 
temperature in all directions. The mean intensity of the 
approx. 1 cm diameter beam as seen by the approx. 
2 mm diameter cloud was measured to be 4 x 10~ 3 I sat 
using a 1 mm diameter pinhole in front of a calibrated 
photodiode unless otherwise noted. The atomic ensem- 
ble is detected by absorption imaging using a 12 Bit CCD 
camera. The resulting optical density profile of the cloud 
is fitted with two orthogonal Gaussian profiles meeting at 
the center of the cloud. The number of atoms is extracted 
from this fit using the peak optical density and the size. 
We perform time-of-flight (TOF) sequences to obtain the 
temperature in axial and radial direction. Densities are 
derived using the size of the cloud in the trap obtained 
from the time-of-flight fit and the number of atoms. The 
error bars shown in the figures are the square root of the 
diagonal elements of the covariance matrix for the fitting 
parameters obtained from a least square fit. No system- 
atic errors have been included, unless otherwise noted. 



IV. EXPERIMENTAL RESULTS 

A. Dynamics 

In this section, we present experimental results on the 
dynamics of the cooling process. In Fig. [| we have plot- 
ted the temperature of the cloud in axial (z) and radial 
(y) direction after a variable cooling time. We start in 
both directions with a temperature of approx. 1 mK. 
In axial (cooling laser beam) direction we see a fast de- 
crease in temperature. But also in the radial direction, 
where no laser beam is applied, cooling is observed on a 
short timescale. The initial decrease in temperature is 
very well fitted by an exp onential decay in both cases, as 
is expected from Eq. |l6|. The cooling time constants de- 
rived from the fit are rf° l ss 10 ms, corresponding to the 
fastest possible cooling time of j^- and Ty°° l ~ 50 ms. In 
fact, we have never observed axial cooling time constants 
significantly lower than 10 ms. The minimum steady 
state temperatures in Fig. |3| are = 334 ± 7 /iK and 
T£° = 124 ± 3 fiK, corresponding to 2.7 and 1 times the 
ID Doppler temperature of To = 124/iK. Due to den- 
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FIG. 3: Typical evolution of the temperatures in axial (z) 
and radial (y) direction as a function of cooling time. The 
intensity of the cooling laser was 4 x I sa t- 
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FIG. 4: Evolution of the number of atoms and the phase space 
density for the same data set as shown in Fig. |S| 

sity dependent heating effects in the light field J25|, p6[ , 
the minimum axial temperature of T™ m ps 90/iK is not 
reached. Fig. |] shows the evolution of the number of 
atoms and the phase space density during the cooling 
in a double logarithmic plot to resolve the fast initial 
dynamics. The number of atoms stays constant within 
experimental uncertainty beyond 300 ms at which time 
the steady state temperature is reached. We gain more 
than a factor of 80 in phase space density and increase at 
the same time the peak density by more than one order 
of magnitude from 0.4 x 10 10 to 5.6 x 10 10 cm~ 3 . 

In the experiment shown in Figs. |3| and ^, we contin- 
ued cooling after reaching steady state. Fig. ^ shows a 
strong reduction in the number of atoms for cooling times 
longer than 500 ms. Atoms are lost from the trap pre- 
sumably by radiative escape or fine structure changing 
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FIG. 5: Radial cooling rate coo i as a function of measured 

single beam cooling light intensity. The corresponding cooling 
rates in axial direction are larger by a factor of at least 5 for 
each data point. 



collisions Ed, 23 in the cooling light. At the same time 



the temperature increases in radial direction, whereas it 
stays constant in axial direction. This is evidence for a 
density dependent cooling effect in radial direction. A 
lower number of atoms results in a reduced optical den- 
sity and, according to Eq. [l4|, to a higher steady state 
temperature T^°. We will come back to this effect in the 
discussion of the steady state results. 

We have already pointed out in section 2B that the 
cooling rate of the radial degrees of freedom via an- 
harmonic mixing or elastic collisions between the atoms 
should be independent of cooling light intensity, provided 
cooling in axial direction is much faster. In Fig. | the re- 
sults of a cooling experiment with different cooling light 
intensities are plotted. For each data point we have per- 
formed an experiment analogous to Fig. || and fitted 
the initial exponential decay of the temperature to ob- 
tain the cooling rate. The corresponding rate constant 
for the axial direction has been verified to be larger by 
at least a factor of five for each data point. One can 
clearly see a linear dependence of the radial cooling rate 
on light intensity, ruling out anharmonic mixing or elastic 
collisions as the major cooling mechanism. We have also 
performed an experiment in which we stopped cooling af- 
ter the axial direction reached steady state, thus produc- 
ing a highly anisotropic temperature distribution. The 
subsequent cross-dimensional relaxation of the tempera- 
tures driven by elastic ground state collisions occurred on 
a much longer timescale than the observed cooling time 
constants when the light was on. 

Increasing the laser intensity beyond the values shown 
in Fig. [| results in a strong distortion of the cloud due 
to light pressure forces and intensity imbalance between 
the two counterpropagating cooling beams. This causes 
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FIG. 6: Radial steady state temperature as a function of cool- 
ing light intensity. The data is fitted by Eq. The data is 
from the same experimental run as in Fig. g. 



parasitic heating and increases the achievable tempera- 
tures. 

From the slope of the linear fit to the data in Fig. || we 
can extract k v according to Eq. [w| Taking into account 
that the average intensity incident on an atom is a factor 
of 3 lower than the measured intensity due to absorption, 
we obtain k v = 0.023 for a detuning of A po i = —0.8. To 
compare this experimentally determined value with the- 
ory, we have integrated Eq. ^ using the steady state size 
of the cloud a z — 700/im and <j y — 410/im. The theo- 
retical result of K y = 0.1 is larger by more than a factor 
4. This deviation supposedly arises from uncertainties in 
the effective intensity as seen by the atoms. This includes 
a systematic error in the measured intensity, uncertainty 
in the averaged intensity due to absorption of photons 
from the laser beam and averaging effects for the detun- 
ing due to magnetic field curvature. 



B. Steady State Temperatures 

In this section, we present experiments that elucidate 
the different effects influencing the steady state temper- 
ature and compare them to theory. 

In the simplified model used to derive the steady state 
temperatures in axial and radial direction given by Eqns. 
|l3| and [l4[ the latter are independent of cooling light 
intensity, since both, the heating rate 1/T^ at and the 
cooling rate 1 /t^°° 1 are linear in intensity. Accounting for 
an additional constant heating rate R in the differential 
equation (^|) , the steady state temperature for a detuning 
of 5 = —0.5 becomes 



y ~ 



(1 



R , 



- D 



(17) 



3000 
2500 
2000 
H500- 
1000 
500 






S data 




theory 


















■ ■ 




ft 









and k* = k/OD v , the steady state temperature reads 
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This equation is plotted as a solid line in Fig. [?] using 
the results for R and K y (normalized to steady state op- 
tical density of 4.5) from the fit to the data in Fig. ||. 
The agreement with the experimental data is surprisingly 
good, since K y and k are not strictly linear in the optical 
density (see Eqns. ||, ||) and the heating rate depends on 
temperature and density. 



CONCLUSION 



FIG. 7: Steady state temperature in radial direction as a 
function of optical density (OD) in this direction. The error 
bar resembles the typical error from the temperature fit to 
the time of flight data. The theory curve is a plot of Eq. |l8| 
using the results of the fit from Fig. pi 



In Fig. U we have plotted the measured steady state tem- 
perature T£° for different cooling light intensities. The ht 
to Eq. using the measured detuning of 6 = —0.8 and 
n = 0.24 calculated from Eq. |, gives n y = 0.11 ± 0.002 
and R = 2.6 x 10~ 26 ± 10% Js _1 . This result is in ex- 
cellent agreement with the integrated value for n y — 0.1 
obtained from Eq. ^. From Fig. [(] we can also estimate 
a minimum radial temperature of T£° w 300/xK which is 
in good agreement with the predicted value using Eq. |lj 
and Fig. | of T a °° = 2.3T D w 290/iK. Thus, the addi- 
tional, intensity independent heating mechanism slightly 
increases the minimum temperature. We have identified 
this mechanism to be a strong dipolar relaxation process 
in chromium ]27| . 

We have already seen in Fig. || that the steady state 
temperature in radial direction depends on the number of 
atoms in the cloud (increase in temperature for a decreas- 
ing number of atoms at long cooling times). To elucidate 
this effect more quantitatively we have performed an ex- 
periment in which we cool down the atoms to steady 
state as described before. In a second "cooling" stage 
we reduce the cooling laser intensity to 1.6 x 10 -3 I S at 
to slow down the loss of atoms. After different cooling 
times, we measure the temperature and the optical den- 
sity in radial direction. The result is plotted in Fig. |?] 
as a function of optical density in radial direction. We 
observe a decrease in temperature with increasing optical 
density. The effective number of reabsorbed radial and 
total cooling photons, n y and k, respectively, are to low- 
est order linearly proportional to the optical density OD y 
(Eq. |^). Introducing two new constants, k* = n y /OD y 



We have presented a general Doppler cooling technique 
which works for optically dense, trapped atomic samples. 
In experiments with magnetically trapped Chromium 
atoms we have achieved an increase in phase space den- 
sity by two orders of magnitude. To explain our experi- 
mental results, in particular the fast and efficient three- 
dimensional cooling in a one-dimensional molasses, we 
have developed a model which takes into account radial 
cooling due to reabsorption of scattered photons. Our 
experimental findings for the intensity-dependence of the 
cooling rate and the steady-state temperature agree well 
with the theoretical considerations. We could exclude 
thermalization effects via elastic collisions or anharmonic 
mixing as the origin of the radial cooling we observe. Op- 
timum radial cooling is observed in spherical trap config- 
urations and the radial cooling efficiency increases with 
the radial optical density. 

Our technique should be applicable to most atomic 
species which can be laser-cooled provided that a mod- 
erate magnetic offset-field (e.g. Bo = 15 G for 87 Rb) is 
applied to prevent depolarization of the sample. Imple- 
mentation of this technique in a BEC experiment could 
simplify the initial laser cooling stages, since Doppler 
cooling in the trap reduces the figure of merit for load- 
ing of the trap to the number of atoms transferred into 
the trap regardless of temperature. Due to the low light 
intensities needed to cool the atoms, our method might 
also be applicable to magnetically trapped atoms in a 
cryogenic environment (28], [2{|. In our experiment, the 
achieved phase space density for chromium already pro- 
vides good starting conditions for evaporation and only 
excessive atom loss and heating due to dipolar relaxation 
|p7f in the consecutive evaporation stage prevented us 
from reaching quantum degeneracy. 
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